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ABSTRACT 


This study Investigates the steady-state and transient 
response of the squeeze film damper bearing. Both the steady- 
state and transient equations for the hydrodynamic bearing forces 
are derived. The steady-state equations are used to determine the 
damper equivalent stiffness and damping coefficients. These co- 
efficients are used to find the damper configuration which will 
provide the optimum support characteristics based on a stability 
analysis of the rotor-bearing system. The effects of end seals 
and cavltated fluid film are Included. 

The transient analysis of rotor-bearing syst ms is performed 

by coupling the damping and rotor equations and integrating 

forward in time. The effects of unbalance, cavitation and 
retainer springs are included in the analysis. 

Methods of determining the stability of a rotor-bearing 
system under the influence of aerodynamic forces and internal 
shaft friction are discussed. Particular emphasis is placed on 
solving the system characteristic frequency equation and stability 
maps produced using this method are presented. 

The study shows that for optimum stability and low force trans- 
mlssablllty the squeeze bearing should operate at an eccentricity 
ratio of c<0.4. 
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1. INTRODUCTION 


Modern turbomachlnes are highly complex systems. Current 
design trends are producing machines that consist of several . 
process stages joined together. The rotors In these machines 
are highly flexible shafts, often mounted in more than two bear- 
ings, that rotate at very high speeds. It is not uncommon to 
see machines that operate above the second critical speed. As 
a result the system dynamics are very complicated. 

One of the major problems encountered in these machines is 
instability produced by aerodynamic forces on impeller wheels, 

friction in the stressed rotor and hydrodynamic forces in the 
bearings. The instability is characterized by large amplitude 

whirl orbits and often results in bearing or total machine failure. 

It is often aggravated by seals and other external forces 

transmitted to the machine. Production losses from failed machines 

are very high and it may take many months to repair or replace the 

failed unit. In addition operator safety is jeopardized when 

machines fall and occasional loss of life occurs. 

An an example of violent self excited whirl motion in a turborotor 

is Fig. 1-1 which represents a centrifugal compressor designed to operate 
at 52,000 RPM. At the design speed the compressor becomes highly unstable, 
the magnitude of the whirl motion being dependent upon the compressor 
speed and loading. The occurrence of this type of whirl motion has 
caused several catastroplc failures with this class of compressors. 
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ROTOR amplitude, M*V 






From the earliest investigations of rotor instability, it has 
been known that the use of flexible, damped supports has an effect 
on Instability and can eliminate it or alter the speed at which 
it occurs. Recent research has produced a large body of knowledge 
on the use of these supports and their effect on instability. 

The squeeze film damper bearing is one type of flexible 
support that is currently being investigated. This study ex- 
amines the squeeze bearing and through computer simulation shows 
its effects on several rotor-bearing systems. The equations for 
the hydrodynamic damper forces are developed in both fixed and 
rotating coordinate systems. The use of two coordinate systems 
allows for both steady-state and transient analysis of damper 
- performance. 

The steady-state behavior of the damper results in the for- 
mulation of damper stiffness and damping coefficients which can 
be used to ^ize the damper configuration. This is accomplished 
by comparing the coefficients with required values obtained from 
a stability analysis of the rotor-bearing system. 

The transient analysis is very useful in determining the 

i 

response to particular forms of external and Internal 
forces as noted previously. Also the effect of damper retainer 
springs and fluid film cavitation can be found. The transient re- 
sponse is found by tracking the Journal motion forward in time 
by integrating the equations of motion under the influence of the 
system forces. 
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2. THEORETICAL ANALYSIS 


2.1 REYNOLDS EQUATION 


The configuration of the squeeze film damper bearing 
is shovm in Figure (2-1) vhere the clearance has been exag- 
gerated. Both fixed and rotating coordinate systems are shovm, 
and the bearing equations are derived for both systems. The 
definitions of the various parameters are listed in the nomen- 
clature section of this paper. 

The basic bearing equation is the Reynolds equation which 
is derived from the Navler-Stokes equations for incompressible 
flow. With the proper bearing parameters the equation for the 
fluid film forces are derived. [1] 

The Reynolds equation for the short, plain journal bearing 
is given In both fixed and rotating coordinates by: 

Fixed coordinates: 


d fh* 3P] _ ... 4. X 3h . 2 ah 

[sir 5z J ■ “j> W * Jt 

Rotating coordinates: 

Iz [ly H] “ “j “ + 2 ~ 


(2-1) 


( 2 - 2 ) 


As shown In Figure (2-1), the angle 6 in the fixed coor- 
dinate expression Is measured from the positive x-axls In the 
direction of rotation whereas the angle 6* in the rotating co- 
ordinate expression is measured from the line of centers In the 
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direction of rotation. The assumptions used In the derivation 
of equations (1) and (2) Include: 

1. The fluid Inertia terms In the Navler-Stokes equations 
have been neglected due to their small magnitude. 

2. Body forces In the fluid film have been neglected. 

• 3. The fluid viscosity Is constant. 

4. The flow In the radial direction has been neglected, 
that Is, the short bearing approximation has been used. 

A comparison of the short bearing solution and the general 
solution of the Reynolds equation solved by a finite difference 
technique for the plain journal bearing under steady state con- 
ditions shows that the short bearing solution is highly accurate for 
a wide range of eccentricities for L/D < 1/4 and Is acceptable for 
L/D values up to 1 If the eccentricity ratio is low. The normal 
design range of the squeeze film bearings will be L/D < 1/2 and 
eccentricity ratios < 0.4. 

2.2 BEARING FORCES IN FIXED COORDINATES ' 

For th^ plain bearing with full end leakage the appropriate 
boundary conditions are: 

(2-3) 

P(e,0) = P(0,L) = 0 

In the fxxed coordinate system the film thickness, h. Is given 
by: 

h « c - X coa 9 - y sin 0 

(2-4) 
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Substituting Into equation (2-1) and Integrating yields: 


P(0 


3u r_9 _ Jl ^ '*’2 I 

' ' f/ L JL ^ 


(2-5) 


The total force components in the x and y directions are found by 

Integrating the pressure over the entire journal 
surface. 


■ 2n L* 




R COS 6 dZ d6 


(2-6) 


2tt l 


- //ne.z 
0 0 


)R sin 0 dZ d0 


(2-7) 


Substituting the expressions for — and ‘ — into the pressure 
equation and integrating around the bearing circumference gives: 



2n 



1 (xsinO-ycos0)-2 (xcos0+ysin0) 

cos 6^ 


f (c - X cos 0 - y sin 0)* 

sin 6j 


P 

The above equation Is applicable to the evaluation of 
the forces developed In the plain Journal bearing as well as 
the squeeze film damper bearing for arbitrary values of journal 
dlspacement, velocity, and shaft and bearing housing angular 
velocities. Hence the analysis can also be used for the general 
floating bush ?jearlng with rotation. 

. For the case of the squeeze film damper where the journal and 
housing ara constrained from rotating, * 0), the force 

expressions become 
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1 

2tt 

• . 

pxl 

_ 1 -2(* cos e + ^ sin e) 

[cos 6 I 

Wj 

2 1 

y (c - X cos e - y sin 0) * 
0 

sin 6 j 


(2-9) 


These non-linear fluid film forces are easily combined with 
the rotor-bearing system dynamical equations providing a complete 
non-llncar dynamical analysis of the system. Because the bear- 
ing force equations are written In fixed Cartesian coordinates 
a transformation from one coordinate system to another Is not 
required. This Is very important for conservation of computation 
time since the bearing pressure profile must be Integrated at 
each time step of the system motion. 


2.3 DAMPER CAVITATION 

If the complete pressure profile Is calculated without re- 
gard to cavitation or rupture of the film, then the bearing pressure 
will be similar to Figure (2-2). This figure represents the three 
dimensional pressure generated in tbe bearing. 

The exact mechanism causing cavitation In fluids Is not 
fully known. It Is known that film rupture is Influenced by gas 
and solid content of the fluid. Recent Investigations have 
shown that a fluid may stand large tensile stresses- [2] , 
and Its ability to withstand rupture Is dependent on Its past 
history. In this investigation It Is assumed that cavitation 
occurs when the pressure In the film drops below ambient pressure. 
The cavltated film then extends over only a section of the bear- 
ing circumference as shown In Figure (2-3) . Recent experimental 
research has shown that cavitation In the squeeze bearing occurs 
In streamers of bubbles which extend around the entire bearing 
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13]. These streamers Initially appear at the center of the bear- 
ing and extend outward as the rotor speed Increases. It Is be- 
yond' the scope of this present paper to analyze this type of 
cavitation effect. Therefore the conventional cavltated film 

Is assumed to occur when P<P where P Is the assumed cavitation 

c c ■ 

pressure. 

When evaluating the Integral of equation (2-9) , negative 
pressures are equated to zero If the film Is assumed to cavitate. 
If the oil supply pressure Is sufficiently high and suitable 
operating conditions exist the film does not cavitate. 


2.4 BEARING FORCES IN ROTATING COORDINATES 


The Reynolds equation In rotating coordinates was given 
by equation (2-2). Assuming steady-state circular synchronous 
precession of the Journal about the bearing center and no axial 
misalignment, equation (2-2) can be Integrated In closed form. 

The resulting equations for the bearing forces give the equivalent 
stiffness and damping of the bearing. 

The force components are: 


ej 


• 





_ -pRL> 

j (4 e sin 6 ' 

+ e cos 0*) 

(cos 

0' 

■ "■ 1 

[ (1 + E COS 

0') ’ 

jsin 

0^ 

/ 


ae' (2-10) 


The limits of Integration, 6] and 0^, define the area over. 


which a positive pressure profile exists and are dependent on 
the type of Journal motion and whether or ‘not cavitation occurs. 

It Is assumed that the damper ig precesslng In steady- 
state circular motion about the origin and therefore e ” 0. 
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The resulting force components are: 


_ -2yRL*cwe_ (2-11) 

^ " C*(l - eM* 

and 

Fq = “PRL^TTeo) (2-12) 

2c’ (1 - 

The force in equation (2-11) appears as a stiffness coefficient 
times a displacement acting in line of the displacement towards 
the bearing center. The equivalent damper stiffness is: 



2yiRL’ eo) 

C* (1 - G*) * 


(2-13) 


Since the journal is processing and not rotating, every point in ' 
the journal has a velocity equal to ew. The force in equation (2-12) 
therefore appears as a damping coefficient times a velocity acting 
in the direction opposite the journal motion. 

The equivalent damper damping is: 


Cq = 

2c* (1 - e*) 3/2 

For the uncavltated film, the components arc given by: 


(2-14) 



(2-15) 


„ _ -pRL*Treo) 

« cMl - 


(2-16) 


It is therefore evident that a complete fluid film does not 
produce an equivalent bearing stiffness but doubles the damping 
of the cavitated film. 
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Although the equations for the damper characteristics were 
derived for a plain damper* they are applicable to other damper 
configurations such as those shown in Figure 2-4. A damper with 
a circumferential oil groove and full end leakage, shown in Figure 
2*'4a, consists of two plain lands of length L/2. The pressure 
drops to atmospheric in the oil groove creating two parallel pres- 
sure profiles corresponding to two parallel dampers of length L/2. 
Thus the maximum pressure in each side of the damper is reduced 
a factor of four and the force on each side is reduced by a fac- 
tor of eight. The net effect of cutting a circumferential oil 
groove in the damper is to decrease the hydrodynamic force and 
both bearing coefficients by a factor of four. 

The damper represented in Figure 2 -4b has a circumferential 
oil groove and end seals to prevent end leakage. Since the pres- 
sure vanishes at the center, the pressure profile is equivalent to 
that for a plain land without an oil groove and without end seals. 
The values for total force, stiffness, and damping are also the 
same. 

As will be shown in later sections, it is often desired to 
design a damper to have stiffness and damping values within a 
certain range. The amount that the clearance can be increased 
may be limited by sealing or other conditions which necessitates 
lowering the damper characteristics by other means. Circumfer- 
ential grooves serve to lower both stiffness and damping of the 
damper without changing the distance between the end seals. The 
characteristics can be easily calculated by summing the effects 
of the lands between the oil grooves. 
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b. Damper with Circumferential Oil Groove 
and Flexible "0" Ring End Seals 


Figure 2—4 Dampers with Circumferential Oil Grooves and 
End Seals 
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The damper equations derived In this section are summarized 
In Table (2-1) . Also Included In the table are the equations f^r 
pure radial squeeze motion. For this type of operation ^ - o and 
results from a purely unidirectional load on the Journal. The 
radial and tangential force components are derived from equation 

O 

(2-10) where only the term containing e In the Integral Is re- 
tained. The pressure equation Is also modified to Include only 

• I 

the c term. The maximum pressure occurs at 6 ■ir for all values 
of Journal eccentricity. Examination of the pressure equation 
reveals that the hydrodynamic pressure Is positive only In 
the region 10' ■ y to These values of 0’ are the limits of 

Integration In equation (2-10) for the cavltated film. 

The table also shows that for purely radial motion no damper 
stiffness Is obtained In either the cavltated or uncavltated 

damper. Thus-lf.- this- type.r>£ motion-exlsts-retalner-sprlngs 

must be Included to provide support flexibility. * 

For the case of circular damper precession, the table 
shows the stiffness and damping of the cavltated film and damp- 
ing of the uncavltated film remain essentially constant for low 
eccentricity ratios. As the eccentricity ratio Increases above 
0.4 there Is a rapid Increase In these properties and they approach 
Infinity as e approaches 1. This variation of stiffness In the 
cavltated film Is very Important. As the eccentricity becomes 
large the support becomes more rigid with a corresponding In- 
crease In the rotor critical speed. If the rotor critical speed 
la Increased above the operating speed, the phase angle between the 
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TYPE OF MOTION 

MAXIMUM PRESSURE 

EQUIVALENT DAMPING 
Ko (Ib/in) 

EQUIVALENT DAMPING 
Co (Ib-sec/in) 

CIRCULAR SYNCHRONOUS 
PRECESSION 

^ ■ ta», £ ■ 0 

-SyL^ojc sin 0 

m 

2c^(l + e cos 0 

El 

where 0 is given by: 

IQ 

w 

- Vi 
2 c 3(1 - c2) 

CAVITATED FILM 



UNCAVITATED FILM 

(1+ECOS0 )cos0 +3e8in^0 -0 
mm m 

0 

’’ ifi ■ 
c3(l - e2) ' 

PURE RADIAL SQUEEZE 
MOTION 

-3pL^e cos 0 

m 

0 

‘ liRL3[ir-cos* (e)l(2c^l) 
c3 (1 - 

♦ - 0, e ft 0 
CAVITATED FILM 

UNCAVITATED FILM 

2c^(l + e cos 0 * 

■ 


9 - ir 

m 

0 

wRL3ff(2c2 + 1) 
c^(l - e^) 


Table 2-1 


Summary of Equivalent Stiffness and Damping Coefficients 
for Squeeze Film Damper Bearings. 








rotor unbalance vector and amplitude vector becomes less than 
90**. Wlien this condition occurs the force transmitted through 
the support structure will always be greater than the unbalance 
load. With an uncavltated film this problem does not occur because 
no damper stiffness is generated. To obtain the stiffness re- 
quired to stabilize a rotor (see Section 3) it is necessary to 
use retainer springs In the support bearings. 

One of the most significant parameters affecting damper per- 
formance Is the length to clearance ratio. The stiffness and 
damping coefficients vary as (L/C)^ an^ therefore either doubling 
the bearing length or decreasing the clearance by 1/2 will In- 
crease the coefficients by a factor of 8. 
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3. ROTOR-BEARING STABILITY AW) PMOTl MULtSIS 


3.1 ROTOR-BEARING STABILITY 

After Jcffcott's [5j analysis In 1919 of the single mass 
flexible rotor on rigid bearings , manufacturers began producing 
light, flexible rotors operating above the first critical speed. 
However, some manufacturers encountered severe operating diffi- 
culties with some of their designs. These machines underwent 
violent whirling while running above the critical speed and often 
felled. 

Experimental and analytical Investigations by Newkirk and 
Kimball [6] [7J revealed that the whirl Instability was not caused 
by unbalance In the rotor, but by Internal shaft effects such as 
Internal friction. Kimball theorized that forces normal to the 
plane of the deflected rotor could be produced by alternating 
stresses In the me cal fibers of the shaft. In light of this theory, 
Newkirk concluded also that the same normal forces could be pro- 
duced by shrink fits on the rotor shaft. By Incorporating these 
forces In Jeffcott's model Newkirk showed that the rotor ceuld be 
unstable above the rotor critical speed. 

Further investigation by Newkirk showed cases of rotor In- 
stability which were not produced by shaft effects but by effects 
In the Journal bearings. [8] One cause of Journal bearing In- 
8ti.!'jllty was later shown to be due to lack of radial stiffness 
in the bearing and the instability occured at twice the rotor 
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critical speed. These Instabilities vere especially common In 
lightly loaded rotors and larger bearing loads tended to promote 
stability. The effect of the larger loads Is to cause cavitation 
of the fluid film which results* In a radial stiffness component 
of the bearing forces being produced. 19], [10], [11] 

In 1965 Alford reported on the effects of aerodynamic forces 
on rotors [12]. He showed that these forces couple the rotor 
equations of motion and can produce Instability. He also noted 
that labyrinth seals and balance pistons also produce forces that 
can promote Instability. 

Recent Investigators Including Gunter, Kirk and Choudhury 
[13], [14], [15] have analyzed' the effects of support flexibility 
.^nd damping on reducing rotor Instability produced by the forces 
just described. As a result they have derived stability criteria 
for determining the necessary support characteristics. 

One of the most general methods for determining rotor 
stability Is to derive the characteristic frequency equation of 
the system. The stability Is given by the roots of this equation. 
The real part of the root corresponds to an exponentially in- 
creasing or decreasing function of time. Thus a positive real 
part Indicates instability whereas a negative real part'* In- 
dicates a stable system. This type of stability analysis of a 
rotor-bearing system therefore requires that the characteristic 
equation be known. This equation Is not always easy to obtain. 

The characteristic equation Is derived from the homogeneous 
second order differential equations of motion of the system [15]. 
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By ns.sumlng solutions of the form 

Xj, “ 1“1, 2, ...»n 

and differentiating, the equations are substituted back Into the 
equations of motion. This produces a matrix known as the character- 
istic matrix. The determinant of this matrix gives the characteristic 
equation, a polynomial of degreee 2n in X, where n Is the number 
of degrees of freedom of the system. 

The computer program SDSTB [16] was used to produce the 
stability maps shown In this section. The program calculates 
the characteristic • equation for a three-mass symmetric flexible 
rotor mounted In journal bearings and supported in squeeze film 
damper bearings. The rotor-bearing model is shown in Figure (3-1). 

The rotor Is assumed to remain stationary In the axial direction 
so the rotor has six degrees of freedom and the characteristic 
equation is therefore of degree twelve. The characteristic matrix 
is shown in Figure (3-2) . The determinant of this matrix gives 
the characteristic equation. The unknown variable In this 
equation is X, the natural frequency of the system. An examina- 
tion of the characteristic matrix shows that the coefficients 
of X are functions of the rotor and bearing properties as well as in- 
ternal shaft friction, absolute rotor damping and aerodynamic cross 
coupling. The natural frequencies and stability of the system 
are found by finding the roots of this equation. 

The Journal and support bearing characteristics can either be 
inserted directly as linear coefficients or they may be calculated 
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Figure 3-2 Characteristic Matrix For Three-Mass Model . 

Including Aerodynamic Cross-Coupling Internal 
Shaft Friction and Absolute Shaft Damping 
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in the program from the bearing parameters by solving for the 
equilibrium positions of the journal and support. These charac- 
teristics are non-llnear functions of the journal eccentricity. 

The stability maps In this chapter were produced with the linear- 
ized journal and support bearing characteristics given as Input 
data to the program. The assumption of linear bearing charac- 
teristics Is useful because for low eccentricity the characteristics 
do not vary greatly with changes In eccentricity. This assumption 
allows a large savings In computer time. If the non-llnear charac- 
teristics are calculated, the amount of computer time Increases 
because an iterative procedure Is used to find the equilibrium 
position. 

As an example of how a stability map Is produced, consider 
the following system: 

ROTOR CHARACTERISTICS 


ROTOR WEIGHT 

W2 

at 

675 

lbs 

JOURNAL WEIGHT 

WJ 

- 

312 

lbs (each) 

SUPPORT WEIGHT 

W1 

m 

15 

lbs (each) 

SHAFT STIFFNESS 

KS 

B 

280000 

Ib/in 

SHAFT DAMPING 

CS 

m 

.10 

Ib-sec/ln 

INTERNAL DAMPING 

Cl 


.0.0 

Ib-sec/ln 

ROTOR SPEED 

N 


10000 

RPM 


BEARING CUARACTIJIISTICS 

K - 1.287 X 10^ Ib/ln 

XX 

K - 1.428 X 10^ Ib/ln 

yy 

C “ 1200 Ib-sec/ln 

XX 
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t; “ 1290 Ib-sec/in 


yy 




K 

. K 

m 

0.0 Ib/in 

xy 

yx 


c • 

- C ' 

a 

0.0 Ib-sec/in 

xy 

yx 




Two values of aerodynamic cross coupling were selected, Q 
20000 Ib/ln and Q - 100,000 Ib/ln. For each value of Q, several 
values of support stiffness were selected ranging from 50,000 
Ib/ln to 500,000 Ib/ln. For each value of support stiffness a 
range of support damping values from 0 to 10000 Ib-sec/ln was 
used. Using this method a stability contour was found for a 
given value of aerodynamic cross coupling and support stiffness. 

The rotor and bearing characteristics remained unchanged. 

Figures (3-3) and (3-4) show the stability maps for the 
above system for the two values of aerodynamic cross coupling. 

There Is an Intermediate range of support damping values for which - 

the sytem ls_stable for a given value ^pf the support. stiffness. 

As the stiffness Is Increased the system becomes less stable. 

With Q “ 2(^000 Ib/ln the optimum amount of damping ranges from 
500 to 2500 Ib-sec/ln as the stiffness Increases from 50000 to 
500000 Ib/ln. For damping less than 100 Ib-sec/ln the system Is 
unstable for all values of stiffness. The same Is true If the 
damping exceeds 10000 Ib-sec/ln. 

For Q " 100^000 Ib/ln the optimum damping Is 1000 Ib-sec/ln 
" and does not shift over the stiffness range selected. When the support 
stiffness reaches 25(^000 Ib/ln the system Is unstable for all values 
of damping. 
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REAL PART- ROOT OF CHARACTERISTIC EQUATION 



SUPPORT DAMPING (Ib-sec/jn) 


Figure 3-3 Stability of a Flexible Rotor With Aerodynamic Cross 
Coupling (Q = 20,000 Ib/ln., N = 10,000 RPM) 
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Figure 3-4 Stability of a Flexible Rotor With Aerodynamic Cross 
Coupling (Q = 100,000 Ib/ln., N = 10,000 RPM) 


3.2 STEADY-STATE ANALYSIS 


The linearized stability maps just discussed provide In- 
formation on the support characteristics needed to promote 
stability in a given rotor-bearing system. There remains the 
problem of relating these characteristics to the actual damper 
system. The squeeze bearing equations derived In section 2 In 
rotating coordinates are used to determine the preliminary 
bearing design. As noted In section 2, these equations were de- 
rived assuming steady-state circular » synchronous precession of 
the journal. 

The damper characteristics, stiffness*, damping and pressure 
are functions of the amplitude of the Journal orbit, fluid vis- 
cosity and damper geometry. The addition of oil supply grooves, 

I 

end seals and cavitation affect the damper characteristics. 

The steady-state equations have been programmed on a 'digital 
computer. This program, SQFDAMP, analyzes three basic bearing 
conf Igur a t Ions : 

1. Plain damper no oil supply groove or end seals. 

2. Damper with oil supply groove but without end seals. 

3. Damper with both oil suprly gr'^ove and end seals. 

Both cavltated and uncavltated flupd films can be analyzed. 

Figures (3-5) - (3-7) show the characteristics for a damper 
being considered for the 675 lb rotor system described earlier. 
The damper has an oil supply groove and end seals, and the fluid 
film Is assumed to be cavltated. The damper parameters are, 
length, 1.0 Inches, radius, 1.2 Inches and fluid viscosity 10 
mlcrorcyns. 
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QTXT c QTXT 2 OTXT . T DTXT 

INI/33S-91J - 03 'ONIdWdQ ’ 


SQUEEZE FILM DAMPER 
OIL SUPPLY GROOVE AND SEALS 
. CAVITATION 


Excessive Damping 




Design's. 
Region 
) - 100,000 


Insufficient Damping 


N= 10000.0 RPN. 

R: 1.20 IN. 

L= 1.00 IN. 

PS= 0.0000 P5I. 
Mj=10. 000 NIOtORCrNS 


.M .9 .6 1.0 

ECCENTRICITY - (DIM) 


Figure 3-*S. Damping Coefficient for Squeeze Film Bearing 

With Cavltated Film - End Seals auu Oil Supply 
Groove Included. 
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STIFFNESS, KO - (LB/INl 
1X10 » 1X10 * 1X10 S 1X10 


SQUEEZE FILM DAMPER 
OIL SUPPLY GROOVE AND SEALS 
CAVITATION 
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MRXIMUM PRESSURE - (PSIl 


SQUEEZE FILM DAMPER 
» OIL SUPPLY GROOVE fiNO SEALS 



mW 

ECCENTRICITY - [DIM) 

Figure . Maximum Pressure for Squeeze Film Bearing with 
Cavltated Film - End Seals and Oil Supply 
Groove Included. 
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For the case where Q ■ 20,000 lb. /In. It was determined that the 
optlimim support damping Is about 500 lb. -sec. /In. , and the support stiff- 
ness should be less than 100,000 lb. /In. Because it Is desirable to 
keep the eccentricity ratio of the damper low. Figures (3-5) and (3-6) 
reveal that this damper will provide the necessary stiffness and damping 
characteristics with a clearance of about A mils at an eccentricity ratio 
of c >= .10 to .20. This corresponds to a damper orbit of O.A to 0.8 mils 
amplitude. The maximum hydrodynamic pressure in the damper is about 100 
psl for this clearance. If the fluid film does not cavitate, the re- 
sulting damping characteristics are doubled. A slightly larger clear- 
ance, 5.0 mils, will produce the optimum damping. However, because the 
uncavitated film does not produce an equivalent stiffness, retainer 
springs must be incorporated in the damper. If the end seals are flexi- 
ble, the required spring rate may be obtained from them. 

One advantage of the uncavitated film is that if the journal 
eccentricity ratio should become very large, there is no rise in stiff- 
ness that could cause the system to become unstable or raise the critical 
speed above the operating speed. 

If Q • 100,000 lb. /in., the optimum damping is 1,000 Ib-sec/in, 
and the effective damper stiffness developed by the combined acting of 
the retainer spring support and the damper hydrodynamic action cannot 
exceed 100,000 lb. /In. In order to achieve the necessary damping 
required, a low clearance damper of the order of 3 to A mils Is needed. 
However, since the total support stiffness cannot exceed 100,000 lb. /in., 
then with a 3 mil clearance, the eccentricity of the damper cannot exceed 
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0.15 roil. With a 3 roll clearance and the damper operating with a 
synchronous orbit of c ■■ 0.23 about the daroper origin, a stiffness 
of 250,000 lb. /In. will be generated. This high value of stiffness 
therefore will prevent the rotor system from achieving stabilization. 

If the aerodynamic cross-coupling on the nachlne Is 100,000 
lb. /In., then the squeeze film damper must be sized very carefully 
for the application because of the large aerodynamic loading on the 
machine. For example, on Figure 3-6 a line is drawn at the stiffness 
value of 100,000 lbs. /in. The damper cannot be operating above this 
level because of the excessive stiffness that would be generated by 
the damper. The Intersection of the stiffness value of 100,000 with 
the clearance curves of 3, 6, and 10 mils produces the maximum allowable 
operating eccentricity for this design of damper. On Figure (3-5), the 
line of 1,500 is drawn on the figure. The operation of the damper 
cannot be above this level as this would produce excessive damping In 
the system. Likewise the line at 350 Is drawn on the figure, and the 
damper cannot operate below this limit as there would be insufficient 
damping produced In order to stabilize the rotor. Thus the damper must 
operate so as to produce a range of damping coefficients between 350 and 
1,500 Ib-sec/in If the rotor is to be stabilized with this high value 
of aerodynamic cross-coupling. Next the points of Intersection from 
Fj.gure 3-6 where the stiffness line of 100,000 lb. /In. Intersects the 
clearance curve Is superimposed on Figure 3~5> and this line Is drawn. 

The operation above this line will produce excessive stiffness In the 
daroper. The region bounded by the excessive stiffness line and the 
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damping lines produces the allowable design region. For this case of 
damper under consideration, it can be seen that there Is only a small 
allowable design region between 2 ^ to A mils bearing clearance with 
a maximum eccentricity of .2. If the damper does not operate in this 
region, then proper stabilization of the rotor configuration will not 
be achieved. If the damper system Is to be designed to stabilize a 
Q value of 100,000 lbs. /in., It Is therefore desirable to study other 
damper lengths and clearance values In order to obtain a damper con- 
figuration which has a larger design region. If the Q value for 
design Is only required to be 20,000 lbs. /in., then the design region 
is greatly extended, and the permissible damper clearance can vary 
from 3 to 6 mils with a maximum eccentricity of 2.6. 
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4. ROTOR TRANSIENT ANALYSIS 


In the previous section, the design characteristics of the damper 
were presented based upon linearized rotor bearing theory. If is often 
Important to evaluate the rotor dynamical behavior with and without a 
squeeze film damper under operating conditions. In order to gain an 
understanding of the type of rotor motion generated in the system under 
the action of aerodynamic cross-coupling or internal rotor friction, 
the rotor, as illustrated in the sample problem, was run with a transient 
program to simulate the motion at the rotor center and the bearing loca- 
tions. 

Figure 4-1 represents the journal motion of a turborotor mounted 
in tilting pad bearings operating at 10,000 RPM. The rotor mass center 
has a small unbalance of .5 mils eccentricity. The rotor is operating 
with an internal friction factor of 20 lb. -sec. /in. The behavior of the 
rotor is similar to the case where Q equals 20,000 lb. /in. If noncon- 
tacting probes were placed at the bearing locations to monitor the rotor 
motion, there would probably be little reason for concern as the total 
rotor amplitude is approximately 1/2 of a mil. A frt.iuency analyzer, 
however, would indicate that there is a fractional frequency whirl motion 
present in the journal orbit. 

However, if noncontacting probes were placed at the rotor center, 
then a large whirl orbit of over 4 mils would be detected as shown in 
Fig. 4 t- 2. Therefore, it is seen that the rotor motion is 8 times the 
magnitude of the bearing motion for this particular configuration. Even 
If the rotor speed increases while the aerodynamic cross-coupling or 
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JOURNAL MOTION 

CYCLES 0.00 THROUGH 10.00 
N . 10000 RPM EU = 0.500 MILS 

W2 « 675.00 LB. MJ « 312.00 LB. 

KS = 280.000 LB/MIL KXX = 1287.000 LB/MIL 

CS = 0.00 LB-SEC/IN KYY = 1428.000 LB/MIL 

Cl = 20.00 LB-5EC/IN CXX = 1200.00 LB-SEC/IN 

Q = 0.00 LB/IN CYY = 1290.00 LB-SEC/IN 

WCX = 3629.54 RPM RCX = 53.49 

WCY 3647.25 RPM flCY = 58.92 

TRDB «= 0.652 < 0.76) 



^l.CQO -0.667 -0.333 -0.000 0.333 0.667 i.OQO 


X-DIR. CMILS3 

Figure 4-1 
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KYY = 

1428.000 LB/MIL 
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20.00 LB-SEC/IN 

CXX = 

1200.00 LB-SEC/IN 

Q = 

0.00 LB/IN 

CYY = 

1290.00 LB-SEC/IN 
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3629.54 RPM 

flCX = 

53.49 

WCT = 

3647.25 RPM 

fiCY = 

58.92 



FU = 

958.473 LB. 



Figure 4-2 
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Internal friction remains constant, then the whirl amplitude will grct-itly 
increase. Sustained rotor operations under these conditions could result 
in the contacting of the rotating clement with the casing with possible 
dire consequences. 

Figure 4-3 represents the rotor system with a squeeze damper support 
system added to It with a design support stiffness of 130,000 lb, /in. and 
a support damping of 150 lbs. -sec. /in. The rotor speed is also taken 
as 10,000 RPM. The rotor system is released, and the centerline rapidly 
spirals out from the origin due to the application of the rotor unbalance. 
After approximately 10 cycles of shaft motion, the initial transient motion 
dies out, and the rotor assumes a stable synchronous whirl orbit with a 
maximum amplitude of 1 mil. Figure 4-4 represents the rotor motion with 
the speed Increased to 16,000 RPM. fhe absolute rotor orbit shown In the 
figure represents the trajectory from 10 to 18 cycles of shaft motion. 

Here again it can be seen that the rotor is approaching a highly stable 
synchronous whirl orbit. In Figure 4-3 and 4-4 it is seen that the rotor 
under the action of low internal friction or aerodynamic cross-coupling 
can be stabilized with a damping value as low as 150 Ib-sec/in. However, 
if the damping vale chosen for the damper is too low, the rotor is unsta- 
ble. Figure 4-5 shows several orbits of the rotor with 50,000 lb. /in. 
stiffness, but only 10 lb. -sec. /in damping. This underdamped system 
goes unstable rapidly. 

Figure 3-4 on the stability of the flexible rotor with the large 
aerodynamic cross-coupling of Q » 100,000 lb. /in. shows that in order to 
stabilize the rotor, the support system must be carefully tuned with a 
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T-DIR. CMILS3 


ABSOLUTE ROTOR MOTION 


N = 

10000 RPM 

EU = 
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N2 = 
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WJ = 

312.00 LB. 
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20.00 LB-SEC/ IN 

CXX = 
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3647.25 RPM 
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58.92 
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Figure 4“ 3 
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RBSOLUTE ROTOR MOTION 
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18.00 
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Figure 4-4 
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NO. 606.7201 


RBSOLUTE ROTOR MOTION 
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ACT = 

61.63 
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Figure 4"5 
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Rupport damping of 1,000 lb. -sec. /In. and a support stiffness of the 
range of 50,000 to 100,000 lb. /In. Figure 4-6 represents the rotor 
motion with an aerodynamic cross-coupling of Q ■ 100,000 lb. /In. and 
the optimum support characteristics of K * 50,000 lbs. /In. and a support 
damping of 1,000 lb. -sec. /In. It can be seen that after the Initial 
transient motion has died out the rotor orbit Is a highly stable synchro- 
nous motion with an amplitude of approximately 1.2 mils. If the damper 
should be overloaded In the vertical direction due to external loads 
and Improper centering, then It Is possible for the damper stiffness to 
Increase greatly In excess of the 50,000 lb. /In. optimum value. Figure 
4-7 represents the situation similar to Figure 4~6 except that the verti- 
cal support stiffness has been Increased from 50,000 to 500,000 lb. /In. 

In this case It can be seen that the rotor motion Is highly unstable, and 
the rotor continues to orbit outward with a large self-excited whirl 
motion. After five cycles of shaft motion, the orbit has grown to over 
7 mils amplitude In the vertical direction. 

The transient rotor and damper motion presented In the previous 
figures was based upon linearized bearing and damper coefficients. In 
the actual damper, the forces generated are highly nonlinear functions 
of the journal displacement and velocity comnonents. In order to take 

this Into consideration, transient orbits of the rotor bearing system 
were run using the complete nonlinear damper forces In the program. It 
has been found that the unbalance level, centering and the magnitude of 
the retainer spring rate have a significant effect on the ability of the 
damper to perform properly. For example. Figure 4-7 represents the 
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NO. 602.7207 


ABSOLUTE ROTOR MOTION 

CYCLES 0.00 THROUGH 10.04 
N = 10500 RPM EU = 0.500 MILS 

H? = 675.00 LB. NJ = 312.00 LB. 

K5 = 200.000 LB/MIL KXX = 1200.000 LB/MIL 
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Cl = 0.00 LB-SEC/ IN CXX = 1200.00 LB-SEC/ IN 
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WCT = 3647.25 RPM RCT = 61.63 

HI = 25.00 LB FU = 1056.716 LB. 

MX = 50.000 LB/MIL CIX = 1000.00 LB-SEC/IN 

KIT = 50.000 LB/MIL CIT = 1000.00 LB-SEC/IN 
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Figure 4-7 
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notion of the squeeze film damper with a retainer spring of 50,000 lb. /in. 
and a clearance of 7 mils. The rotor has an unbalance eccentricity of 
25% of the bearing clearance of 1.75 mils, which creats a rotating load 
of 3,700 lbs. Although the initial starting position assumed was con- 
siderably removed from the final steady-state orbit. It can be seen that 
the initial transient quickly dies out and the damper assumes a synchro- 
nous circular orbit with an eccentricity of radius of .25. The maximum 
dynamic transmisslbllity encountered was .74 which occurred during the 
first cycle of shaft motion. After the steady-state motion is achieved, 
the final dynamic transmisslbllity is greatly reduced and is of the order 
of .37. Therefore it can be seen that the squeeze film damper, in addi- 
tion to stabilizing the rotor motion also greatly assists in the attenu- 
ation of the forces generated by rotor unbalance. 

If the squeeze damper were perfectly linear, a doubling of the 
unbalance force would cause twice the force transmitted to the support 
system. However, the dynamic transmisslbllity of .37 would remain the 
same. In Figure 4-9, the unbalance eccentricity is increased from .25 
to .5. Here it is seen that the dynamic transmisslbilty has Increased 
dramatically from .74 to 2.15. This implies that more force is trans- 
mitted through the damper support system than there would be if the 
damper were perfectly rigid. This represents a situation where putting 
in a damper is worse than putting in no damper at all. Therefore it 
is seen that if the rotor is operating with high unbalance as well as aero- 
dynamic cross-coupling, the nonlinear forces generated in the squeeze 
film damper would be such that the high equivalent damper stiffness 
generated would not permit the damper support system to stabilize the 
rotor. 
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Figure 4-^8 Vertical Unbalanced Rotor in Squeeze Film Bear- 
ing - Effect of Unbalance Magnitude - Unbalance 
Ecdentricity ■ 1.75 Mils ■ 
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Figure 4- 9 Vertical Unbalanced Rotor In Squeeze Flln Bear- 
ing - Effect of Unbalance Magnitude - Unbalance 
Eccentficlty - 3.30 Mils. 
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5. EXPERIMENTAL ROTOR MOTION 


The squeeze film damper bearing has been applied successfully to 
several centrifugal compressors In order to stabilize them from self- 
excited whirl motion. However In each case, the squeeze film damper 
bearing had to be carefully sized for the particular machine In con- 
sideration. The design of the squeeze film damper for one machine may 
not necessarily produce satisfactory results In another. Figure 5-1 
represents a centlfugal compressor before and after stabilization with 
a squeeze film damper. The trace at the left represents the original 
crblt of the rotor operating at a compressor discharge pressure of 
175 psl. The rotor became highly unstable above this discharge of 
pressure. After the squeeze film damper was employed, the full com- 
pressor discharge pressure could be achieved. The resulting whirl 
orbit as shown In the right hand upper figure shows that the rotor 
system Is highly stable with only a small component of self excited 
whirl motion existing In the rotor. 

Figure 5-2 represents the frequency spectrum for various operating 
speeds of a centrifugal compressor In tilting pad bearings which exhibited 
whirl instability at the design operating speed. The rotor was run 
from 0 to 14,000 RPM, and the motion was recorded on tape and analyzed 
through a real time analyzer. At the various speeds, an analysis was 
made of the frequency components of the rotor motion. For example, the 
45^ line drawn on the plot represents the synchronous rotor motion. 

There Is also a small component of twice the operating speed as seen on 
the chart. At a speed of 10,500 RPM a small component of sclf-excltcd 
whirl Instability was detected. The frequency of this component corre- 
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N = ll,300 RPM, FIRST CRITICAL N, =4,300 RPM 



UNSTABLE ROTOR ORBIT 
DISCHARGE PRESSURE 
Pq =175 psig 
LARGE N, COMPONENT 


STABILIZED ROlOR ORBIT 
Pq =650 psig 
SMALL N, COMPONENT 



STABILIZED ROTOR ORBIT STABILIZED ORBIT 

X 2 AMPLIFICATION FILTERED SIGNAL 

UNBALANCE ORBIT ONLY 
N = 11,300 


Figure 5-1. Rotor Orbits of a Turbo Compressor Before and After 
Stabilization by Damper Supports 
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spends approximately to the first critical speed at 4,200 RPM. As the 
speed of the rotor is Increased, the subsynchronous whirl motion grows. 
Operation under these conditions caused a considerable wearing of the 
seals and bearings which required periodic replacement of the components. 

Squeeze film dampers were designed for this compressor using the 
methods developed in the previous sections. Figure 5-3 represents the 

frequency spectrum of the compressor with the squeeze film dampers 
installed. Note that the rotor system is now highly stable and that 
there is no indication of self-excited whirl Instability. The synchronous 
motion at design speed was also reduced considerably by means of trim 
balancing a coupling. It should be noted that there is no way that bal- 
ancing the rotor alone could reduce or remove the self-excited nonsynchro- 
nous component in the rotor. It may also be of Interest to note that the 
two times component has not been reduced by the squeeze film dampers or 
by balancing the rotor. It may well be that the two times component is 

caused by a misaligned coupling, and that a hot aligrauent should be per- 
formed on the machine. 

Oscilloscope traces of the rotor orbits before and after stabili- 
zation are shown in Figure 5-4. Figure 5-4a shows that the total motion 
Including both synchronous and nonsynchronous components of vibration at 
the operating speed of 13,500 RPM. Following stabilization of the rotor 
with dampers. Figure 5-4b shows nonsynchronous motion on the upper trace 
and total motion on the lower trace. The resulting total motion is 
nearly all synchronous because the nonsynchronous component is very 
small. 
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FREQUENCY SPECTRUM SPEED VS SPEED 
WITH SQUEEZE FILM DAMPER - STABI LI Z E D SYSTEM 


N= 13,500 DESIG N SPEED 
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Figure 5-3 





ROTOR ORBITS OF A CENTRIFUGAL 
COMPRESSOR BEFORE AND AFTER STABILIZATION 

N = 13,500 RPM 
N,= 4,300 RPM 



A. BEFORE STABILIZATION-TOTAL MOTION 
LARGE SYNCHRONOUS AND 
NONSYNCHRONOUS WHIRL MOTION 
SCALE* I MIL/MAJOR DIV. 



B. AFTER STABILIZATION WITH DAMPER 

UPPER TRACE- NONSYNCHRONOUS MOTION 
LOWER TRACE - TOTAL MOTION 


Figure 5-4 
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The history of this compressor over a period of years showed slowly 
Increasing vibration levels following each replacement of seals and 
bearings until wear became excessive. No large vibrations, either 
synchronous or nonsynchronous , were observed in the rotor operation over 
an eight month period after the dampers were installed. 
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6. Conclusions 


1. Self excited whirl notion nay be created in turbomnchlncry by 
one or nore of the following effects; aerodynanic cross-coupling 
Internal friction, fluid film bearings and seals, balance pistons, 
and labyrinths. 

2. The use of tilting pad bearing does not guarantee stable operation 
because of other effects as mentioned above. 

3.. Long multi-stage turborotors operating at several times the rotor 
first critical speed may be very susceptable to self-excited whirl 
instability. 

4. With short rigid rotors, stability nay be Improved by a reduction 
of clearance of the tilting pad bearings, whereas with long flexi- 
ble rotors, a reduction of bearing clearance can lead to cata- 
stroplc failure. 

5. Stability may be improved by the following 

a. Reduction In operating speed or power level 

b. Reduction of Instability mechanisms 

c. Increase in effective shaft stiffness 

d. Increase in effective rotor nodal damping 

6. Properly designed squeeze film dampers can Improve stability and 
unbalance response characteristics of turbomachinery by increasing 
the modal damping. 

7. The damper characteristics must be sized for a particular rotor- 
bearing system. 

8. For a given damper design there Is a limit to the level of self- 
excltatlon that It can stabilize. 
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9. The higher the level of aerodynamic cross-coupling the more 
carefully the damper must be tuned to the rotor, 

10, The self excited motion at the center of the rotor may be 
considerably greater than the motion monitored at or near the 
bearings 

11, Exces' ve squeeze film stiffness or damping will reduce the 
effectiveness of the damper, 

12, The squeeze film damper characteristics are highly nonlinear 
functions of the eccentricity ratio and hence the damper may 
not function properly under excessive loading. 
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